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Topological entanglements between polymer chains are achieved in associating protein 
hydrogels through the synthesis of high molecular weight proteins via chain extension 
using a variety of thiol coupling chemistries, including disulfide formation, thiol-maleimide, 
thiol-bromomaleimide and thiol-ene. Coupling of cysteines via disulfide formation results 
in the most pronounced entanglement effect in hydrogels, while other chemistries 
provide versatile means of changing the extent of entanglement, achieving faster chain 
extension, and providing a facile method of controlling the network hierarchy and 
incorporating stimuli responsivities. The addition of trifunctional coupling agents causes 
incomplete crosslinking and introduces branching architecture to the protein molecules. 
The high-frequency plateau modulus and the entanglement plateau modulus can be 
tuned by changing the ratio of difunctional chain extender to the trifunctional branching 
unit. Therefore, these chain extension reactions show promise in delicately controlling 
the relaxation and mechanical properties of engineered protein hydrogels in ways that 
complement their design through genetic engineering. 
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INTRODUCTION 

Control of advanced mechanical properties of hydrogels is cen- 
tral to their applications. In many situations, such as engineering 
simulant materials of articular cartilage or blood vessels (Drury 
and Mooney, 2003), hydrogels are required to support mechanical 
load and maintain structural integrity. Multiple approaches have 
been developed to address this issue, with emphasis on increas- 
ing the network toughness and elastic moduli, such as prepar- 
ing double-networks (Gong et al., 2003; Sun et al., 2012) and 
fabricating composite materials (Haraguchi and Takehisa, 2002; 
Wang et al., 2010a). In some scenarios where dynamic properties 
such as injectability, self-healing, shape memory, and controlled 
degradation/reinforcement are desirable, non-covalent crosslinks 
and stimuli-responsive triggers can be incorporated into gels to 
fulfiU the application requirements (Kloxin et al, 2009; Holten- 
Andersen et al., 2011; Guvendiren et al, 2012; Classman et al., 
2013). However, many methods developed in synthetic poly- 
mer systems cannot be readily applied to protein hydrogels, and 
manipulating the mechanical response of the gels presents a new 
challenge. Because the function of a protein is related to its hier- 
archical structure and the diversity of amino acid functional 
groups, site-specific and orthogonal reactions are often required 
to preserve protein properties in the final gel. In addition, the 
modification reactions can usually only be performed in aqueous 
buffers due to the poor solubility of proteins in organic solvents. 

Currently there are two main strategies for chemically modify- 
ing proteins to manipulate the mechanical properties of protein 
hydrogels. The first strategy is covalently crosslinking protein 
polymer chains by adding crosslinking reagents (Trabbic-Carlson 
et al, 2003; Li et al, 2011; Chung et al, 2012), catalysts (Shen 
et al, 2005; Lv et al., 2010) and/or enzymes (Davis et al, 2010). As 



this method only modifies the amino acid residues participating 
in junction formation, the protein strands can retain most of their 
function (e.g., elasticity and stimuli responsivity). In addition, 
many mechanical properties at equilibrium, such as modulus and 
maximum swelling ratio, can be controlled by varying the dosage 
of the crosslinking agents according to well-known laws of hydro- 
gel physics and network theories (Graessley, 2004; Tanaka, 2011; 
Kim et al., 2013). In a typical crosslinking reaction, however, the 
formation of network imperfections, such as dangling chains and 
inelastic loops, is usually uncontrollable and difficult to quan- 
tify (Zhou et al., 2012). The second strategy is coupling synthetic 
polymers and proteins to prepare chimeric copolymers (ling et al., 
2008; Sahin and Kiick, 2009; Wu et al, 201 1; Classman and Olsen, 
2013; Classman et al., 2013). Interactions leading to self-assembly 
of the polymer or protein component then create nanostruc- 
ture within the hydrogel, which in some cases can be triggered 
by external stimuli. This method allows delicate manipulation of 
hydrogel mechanics, such as erosion rate, toughness and elasticity. 
However, preparing hybrid hydrogels adds challenges in materials 
synthesis and purification. 

Recently, we developed a facile method for introducing entan- 
glements into protein hydrogels through simple chain extension 
reactions, and we exploited this entanglement effect to engineer 
new mechanical responses into the materials (Tang et al., 2014). 
This strategy minimally modifies the protein molecules, only 
extending protein chains by establishing disulfide end linkages, 
yet it creates a drastic enhancement in many mechanical prop- 
erties, including the low-frequency modulus, resistance to creep, 
extensibility and toughness. The disulfide linkages are redox 
responsive, which provides opportunities to regulate network 
mechanics by redox stimuli. 
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In addition to disulfide coupling, several other thiol 
chemistries are promising candidates for coupling high molar 
mass proteins: thiol-maleimide, thiol-bromomaleimide and 
thiol-ene coupling. Thiol-maleimide chemistry enables site- 
specific modification of cysteine residues and has been a popular 
route for constructing protein-based bioconjugates (Canalle et al., 
2010; Stephanopoulos and Francis, 2011). Recently, the thiol- 
maleimide addition has been found to be reversible, offering 
new opportunities to control the degradation of the conjugates 
(Baldwin and Kiick, 2011). Thiol-dibromomaleimide conjuga- 
tion has emerged as another important "click" chemistry (Smith 
et al., 2010; Jones et al, 201 1; Robin et al, 2013), but only recently 
has its potential in the chain extension reaction been appreciated 
(Cui et al, 2013). The thiomaleimide adduct is redox respon- 
sive: the addition of reducing reagents such as P-mercaptoethanol 
(BME) reverses the reaction and recovers the unmodified pro- 
teins. Compared to other coupling chemistries, the thiol-ene 
reaction is the most rapid coupling method, with extremely high 
bimolecular rate constants in the thiyl-alkene addition, rang- 
ing from 10^ to 10^ M"' s"' (Northrop and Coffey, 2012), 
which is 4-5 orders of magnitude larger than thiol-maleimide 
addition (Schelte et al., 1999). Despite potential side reactions 
(Schoneich, 2008), there exist many successful examples demon- 
strating the use of thiol-ene chemistry in direct protein modifica- 
tion (Dondoni et al, 2009; Weinrich et al., 2010; Valkevich et al, 
2012). 

In this work, thiol-maleimide, thiol-bromomaleimide and 
thiol-ene coupling chemistries for application in protein chain 
extension reactions to produce entangled hydrogels are com- 
pared. In addition, the ability to use chain extension points for 
the modification of proteins with poly(N'-isopropylacrylamide) 
(PNIPAM) side chains is established, demonstrating thermally 
responsive mechanical behavior. By using trifunctional chain 
coupling agents, branched proteins are also prepared, and 
the effect of the ratio of di to trifunctional chain coupling 
agent on entanglements is assessed. These experiments demon- 
strate sophisticated control of molecular structure and network 
mechanics in engineered protein hydrogels, yielding the abil- 
ity to control chain topology, chain entanglement, and chem- 
ical functionalization aU through thiol-based chain coupling 
chemistries. 

MATERIALS AND METHODS 
MATERIALS 

Bismaleimide diethylene glycol (la) was purchased from Thermo 
Fischer. Maleimide-PEG(lk)-maleimide (lb) was purchased 
from creative PEGWorks. The average molecular weight of 
the PEG portion determined by NMR was 1082 Da. 
P-cyclodextrin (pCD) and 2,2'-Azobis[2-(2-imidazolin-2-yl) 
propane] dihydrochloride (VA-044) were purchased from Wako 
USA. All other chemical reagents were purchased from com- 
mercial sources (Sigma-Aldrich and VWR) and used as received 
unless otherwise noted. 

CHARACTERIZATION 

'H NMR spectra were recorded in CDCI3 or DMSO-dg using 
a Varian Mercury 300 MHz Spectrometer in Department of 



Chemistry Instrumentation Facilities at MIT. High-resolution 
mass spectrometry (HRMS) data was obtained on a Bruker 
Daltonics APEXIV 4.7 Tesla Fourier Transform Ion Cyclotron 
Resonance Mass Spectrometer. Matrix assisted laser desorption 
ionization mass spectrometry (MALDI) data was obtain on a 
Bruker Omniflex MALDI-TOF Mass Spectrometer. Q;-Cyano-4- 
hydroxycinnamic acid (CHCA) was used as matrix. Gel perme- 
ation chromatography (GPC) was performed on an Agilent 1260 
system equipped with a Wyatt Optilab T-rEX refractive index 
(RI) detector and a Wyatt Mini-DAWN multi-angle Hght scat- 
tering (LS) detector. The mobile phase was DMF supplemented 
with 0.02 M LiBr and the instrument was operated at 1.0 mL/min 
at 70° C. 

PROTEIN Cys-P4-Cys EXPRESSION AND PURIFICATION 

Protein expression and purification by ammonium sulfate purifi- 
cation have been described previously (Tang et al., 2014). In this 
study, the proteins were additionally purified by anion exchange 
chromatography using a HiTrap Q Sepharose HP 5 mL column 
(GE healthcare, WI), eluting with a gradient of 0-500 mM NaCl 
in 6M urea and 20 mM Tris (pH 8.0). A typical isolation yield 
is 120 mg per liter culture and the protein purity was deter- 
mined to be >97% by sodium dodecyl sulfate-polyacrylamide gel 
electrophoresis (SDS-PAGE). 

CHEMICAL COMPOUND SYNTHESIS 
Dibromomaleimide-alkyne 

The synthesis procedure reported elsewhere (Jones et al, 2011) 
was slightly modified. Potassium carbonate (0.89 g, 6.50 mmol) 
was suspended in 20 mL acetone and 2,3-dibromomaleimide 
(1.5 g, 5.90 mmol) was added to the slurry in one portion and the 
reaction was left to stir at room temperature for 5 min. Propargyl 
bromide (80% in toluene, 0.72 mL, 6.50 mmol) was added drop- 
wise to the mixture over 10 min. After 24 h, solvent was removed 
under vacuum and the mixture was redissolved in DCM. Salts 
were filtered and the residue was loaded onto a silica gel column 
vacuum. The crude product was purified by flash chromatogra- 
phy, eluted with 0-2% MeOH in DCM (TLC Rf = 0.78, stained 
with KMn04 solution) to afford 492 mg 4 as a white powder 
(yield 28.5%). 

iH NMR (300 MHz, CDCI3) S 4.38 (d, / = 2.4 Hz, 2H), 4.83 
(t, / = 2.4 Hz, IH). DART HRMS (m/z) calcd for C7H4Br2N02 
[M-|-H] + : 293.8586; found 293.8585. 

Dibromomaleimide-p -CD 

To a 25 mL Schlenk tube was added 4 (58.5 mg, 0.20 mmol), 
mono-6-deoxy-6-azido-P-cyclodextrin (Petter et al., 1990) 
(116.0mg, O.lOmmol), Cu(I)Br (14.3 mg, O.lOmmol) and 5mL 
DMF. The mixture was degassed through 3 freeze-pump-thaw 
cycles. 2,2'-Bipyridine (15.6 mg, O.lOmmol) was added to the 
frozen mixture and the mixture was degassed one more time. 
The click reaction was performed at 30°C for 24 h. Catalyst was 
removed by passing through a short alumina column and the 
product was obtained by precipitation in acetone twice to get 
63 mg 2 as a yellow powder (yield 43.4%). 

iH NMR (300 MHz, DMSO-dg) S 8.06 (s, IH), 5.95-5.60 
(m, 14H), 5.10-4.67 (m, 9H), 4.65-4.42 (m, 6H), 3.79-3.48 
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(m, 28H), 3.43-3.21 (m, overlaps with HOD). MALDI-TOF MS 
(m/z) calcd for C49H72Br2N403Na [M+Na]+: 1473.22; found 
1473.22. 

Bisallyl tetraethylene glycol 

Tetraethylene glycol (5.0 g, 25.7 mmol) was dissolved in 30 mL 
anhydrous DMF. Sodium hydride (3.5 g, 60% in mineral oil, 
87.5 mmol) was added to the mixture in one portion. The mix- 
ture was cooled with a water bath and allowed to stir for 30 min 
at room temperature. After the reaction stopped bubbling, allyl 
bromide (9.5 g, 78.4 mmol) was added dropwise. After 14 h, 
excess NaH was quenched by adding 10 mL saturated NH4CI 
aqueous solution. The mixture was diluted with EtOAc, washed 
with DI water and brine, and dried over MgS04. Then sol- 
vent was removed under reduced pressure. The crude product 
was purified by silica gel column chromatography, eluted with 
50% EtOAc in hexanes (TLC Rf = 0.41, stained with KMn04 
solution) to afford 3.65 g 3 as a slightly yellow liquid (yield 
51.8%). 

'H NMR (300MHz, CDCI3) 8 5.91 (d x d x t, / = 17.1 Hz, 
/= 10.5 Hz, V= 5.7 Hz, 2H), 5.27 (d x d x t, / = 17.1 Hz, 
/= 1.5Hz, /= 1.5 Hz, 2H), 5.18 (d x d x t, / = 17.1 Hz, / = 
1.8Hz,/ = 1.2 Hz, 2H), 4.02 (d x t,/ = 5.7Hz,/ = 1.5Hz,4H), 
3.68-3.58 (m, 16H). DART HRMS (m/z) calcd for C14H27O5 
[M4-H] + : 275.1853; found 275.1851. 

EMP-Adamantane lEMP-Ad) 

The synthesis of 2-ethylsulfanylthiocarbonylsulfanyl-2-methyl- 
propionic acid (EMP) was performed as previously reported (Lai 
et al, 2002; Convertine et al, 2006). To a 25 mL round bottom 
flask was added EMP (179.5 mg, 0.80 mmol), 1-adamantane 
methanol (159.6 mg, 0.96 mmol), 4-dimethylaminopyridine 
(DMAP) (19.5 mg, 0.16 mmol) and 5mL DCM. After all 
reagents were dissolved, N,N'-dicyclohexylcarbodiimide (DCC) 
(247.6 mg, 1.20 mmol) was added in one portion. The reaction 
mixture was stirred at room temperature overnight. The precipi- 
tate was filtered and the solvent was removed under vacuum. The 
crude product was purified by silica gel column chromatography, 
eluted with 30% DCM in hexanes (TLC Rf = 0.30, yellow or 
stained with KMn04 solution) to afford 210 mg 5 as a bright 
yellow solid (yield 70.5%). 

'H NMR (300MHz, CDCI3) S 3.60 (s, 2H), 3.25 (q, / = 

7.5 Hz, 2H), 1.93 (s, 3H), 1.71-1.60 (m, 12H), 1.49 (d, / = 
2.4 Hz, 6H), 1.28 (d, / = 7.5 Hz, 3H). ESI HRMS (m/z) calcd for 
C18H28O2S3 [M-|-H]+: 373.1343; found 373.1324. 

PNIPAM-Ad 

JV-isopropyl acrylamide (NIPAM) was freshly purified by 
sublimation and azobisisobutyronitrile (AIBN) was recrystal- 
lized twice from ethanol. In polymerization, NIPAM (1.09g, 

9.6 mmol), EMP-Ad (29.8 mg, 0.080 mmol), AIBN (2.63 mg, 
0.016 mmol) and 4.8 mL acetonitrile were added to a 25 mL 
Schlenk tube. The reaction mixture was subjected to 3 cycles of 
freeze-pump-thaw to degas oxygen. The reaction was heated at 
60° C for 4h, after which polymers were recovered by precipita- 
tion in diethyl ether. The molar mass of the obtained polymer was 
7.8kg/mol (from GPC, dispersity D = 1.05), 7.4kg/mol (from 



^H NMR endgroup analysis) and 7.0kg/mol (from MALDI-TOF 
MS), respectively. 

EXTENDING PROTEIN CHAINS USING DIFFERENT CYSTEINE COUPLING 
CHEMISTRIES (INCLUDING BRANCHING REACTIONS) 

Thiol-maleimide 

Protein Cys-P4-Cys was dissolved in denaturing buffer (8M urea 
and 100 mM phosphate, pH 8.0) to reach a concentration of 10% 
(w/v). Tris(2-carboxyethyl)phosphine (TCEP) (20 eq.) was added 
to the solution and the pH was adjusted to 7.5. Bismaleimide 
la or lb was dissolved in DMF, and was added to the solution 
(1 eq.). The reaction was stirred at room temperature for 3 days, 
and it was then dialyzed against MilliQ water and lyophilized. The 
long reaction time was chosen based on experimental results to 
maximize endgroup conversion and entanglement in gels. 

Thiol-bromomaleimide 

Protein Cys-P4-Cys was dissolved in denaturing buffer (8 M urea 
and 100 mM phosphate, pH 8.0) to reach a concentration of 10% 
(w/v). TCEP (20 eq.) was added to the solution and the pH was 
adjusted to 6.2. pCD functionalized dibromomaleimide 2 was 
dissolved in DMSO, and was added to the solution (1 eq.). The 
reaction was stirred at 4°C for 7 days, and it was then dialyzed 
against MilliQ water and lyophilized. 

Thiol-ene 

Protein Cys-P4-Cys was dissolved in denaturing buffer (8 M urea 
and 100 mM phosphate, pH 8.0) to reach a concentration of 10% 
(w/v). TCEP (20 eq.) was added to the solution. The mixture was 
left to stir at room temperature for 4 h, and it was then dialyzed 
against MilliQ water and lyophilized. After reduction, >99% of 
the proteins were in the monomeric state, assessed by SDS-PAGE 
(Supplementary Figure 2). Reduced proteins were hydrated in 
100 mM sodium phosphate buffer (pH 7.6), and mixed with 3 
(1 eq.) and VA-044 (0.2 eq.). The final concentration was adjusted 
to 20% (w/v). Hydrogel samples were loaded on the rheometer 
after hydration for 2 days. The reaction was triggered by heating 
at 60°C for 3 h. The moduli were monitored at 25°C to ensure 
that the steady state was reached before measurement. 

In branching reactions, the triftmctional crosslinker 2,4,6- 
triallyloxy-l,3,5-triazine 6 was used in place of the bifunctional 
crosslinker 3. The total number of alkene groups was kept in 1:1 
molar ratio to thiol groups, while the amount of 6 was varied from 
0-100% of the total mole fraction of alkene-containing oligomer. 

RHEOLOGY 

Rheology experiments were performed on Anton Paar MCR 
301, 501, and 702 rheometers using a cone and plate geome- 
try (25 mm diameter, 1° cone and 50 fim truncation gap) or a 
parallel plate geometry (10 mm diameter and 300 /xm gap, only 
for measurement with pCD-functionalized protein hydrogels). 
The quantitative measurement of hydrogel mechanics was not 
affected by the choice of rheometers. Lyophilized proteins were 
hydrated in 100 mM sodium phosphate buffer (pH 7.6) to a 
final concentration of 20% (w/v). Hydrogels were kept at 4°C 
for 2 days to allow complete hydration. In order to minimize 
dehydration during measurement, the edges of hydrogel samples 
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were coated with mineral oil. After loading onto the rheometer, 
samples were heated from 25°C to 90°C and cooled to 25°C at 
5°C/min. Unfolding of the coiled-coil domains at high temper- 
atures allowed rapid stress relaxation within the gel to eliminate 
any shear history, but no sol-gel transition was observed over this 
temperature range (Supplementary Figure 1). Frequency sweep 
measurements were performed at 1% strain in the linear vis- 
coelastic regime (LVE). In creep experiments, a 25 Pa load was 
exerted on hydrogels for 2 h, and the load was removed to monitor 
the recovery behavior for 2 h. The creep-recovery experiment was 
repeated at 50 Pa load to ensure the deformation was independent 
of the applied load. 

RESULTS AND DISCUSSION 

COMPARISON OF VARIOUS CHEMISTRIES ON THE CHAIN EXTENSION 
REACTION 

Chain extension of cysteine end-capped proteins can be achieved 
by applying thiol-maleimide, thiol-bromomaleimide and thiol- 
ene coupling chemistries under appropriate reaction conditions 



(Scheme 1 and Figure 1). All chemistries require reduction of 
existing disulfides to recover reactive cysteine residues. In order to 
achieve significant changes in mechanical properties, the conver- 
sion of these macromolecular polycondensation reactions must 
be relatively high to produce proteins of high molecular weights. 
In practice, this is achieved by controlling stoichiometry of the 
reacting species. A theoretical full conversion of endgroups can 
only be obtained when the ratio of cysteines to alkenes is 1:1. 
Running the reaction for extended periods is also useful to reach 
high conversions. 

Each chemistry offers its own distinct advantages and disad- 
vantages for chain extension. In the thiol-maleimide reaction, 
pH control is critical to minimize amine-maleimide coupling 
while achieving fast conjugation (Hermanson, 2008). In addition, 
dimaleimide reagents are commercially available with varying dis- 
tance between the two maleimide groups, and this spacer length 
may affect the polycondensation reaction due to their subtle dif- 
ferences in solubility or by altering the propensity to form bridges 
or primary loops (Button et al., 1994). To study the effect of the 
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Scheme 1 | Cysteine coupling chemistries to extend protein chains. (I) Thiol-maleimide conjugation; (II) thiol-dibromomaleimide conjugation; and 

(III) thiol-ene chick chemistry. 
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FIGURE 1 1 Analysis of the chain extension products by SDS-PAGE. 
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distance between two conjugation sites, bismaleimide with two 
different oligo ethylene glycol spacer lengths la and lb are used 
to conduct chain extension. It is found that the spacer length does 
not greatly affect the chain extension when the coupling reagents 
can be solubUized and homogeneously dispersed in the reaction. 

As for thiol-dibromomaleimide coupling, controlling the 
hydrolysis of the maleimide ring, especially the monothio adduct, 
is crucial to attain significant chain extension. If hydrolysis 
happens on the dithio adduct, transformation from maleimide 
to maleamic acid does not affect chain extension, and it 
only results in loss of the reversibility of the thiomaleimide 
adduct (Supplementary Scheme 1). In contrast, hydrolysis of the 
monothio adduct can potentially limit the conversion, as the 
number of reactive bromomaleimide functional group decreases 
(Supplementary Scheme 1). Maleimide hydrolysis can be regu- 
lated in many ways, including changing temperature, pH (Ryan 
et al, 201 1) and the electron density distribution in the maleimide 
structure (Nathani et al., 2013). In our hands, the first two param- 
eters are optimized: a low temperature (4°C) and a slightly acidic 
buffer condition (pH 6.2) are chosen to be the reaction condi- 
tion, under which moderate chain extension is achieved. After 
chain extension, the protein is exposed to 1000-fold excess BME, 
and ~ 98% of the proteins are converted to the monomeric 
state (Figure 1), which demonstrates the reversibility of the 
dithiomaleimide adduct. 

In the thiol-ene coupling strategy, bis allyl compound 3 is 
chosen as the chain extender to prevent homopolymerization 
of alkenes that can occur when (meth)acrylate groups are used 
(Hoyle and Bowman, 2010), and water soluble VA-044 is selected 
as the thermal initiator for its low decomposition temperature. As 
the thiol-ene reaction is tolerant of oxygen (Hoyle et al., 2010), 
no cumbersome degassing procedure is required, enabling its 
convenient use in applications. Upon heating at 60°C for 3 h, 
a large fraction of high-molecular-weight proteins is formed. 
While thermally initiated chain extension is used here in order to 
ensure homogeneous reaction through a concentrated solution, 
it is also possible to perform photoinitiated chain extension using 
this chemistry, especially in occasions when spatial and temporal 



control of gel mechanics is required. TCEP needs to be removed 
prior to the thiol-ene reaction; otherwise the degree of chain 
extension is fairly low (Supplementary Figure 2). The deleterious 
effect of TCEP may be due to the desulfurization of cysteines cat- 
alyzed by TCEP during free-radical-based reactions, an effect that 
has been observed under similar reaction conditions (Wan and 
Danishefsky, 2007). 

Cysteine coupling chemistries exhibit different performances 
in the chain extension reaction, and the differences in the molec- 
ular weight distribution of the chain-extended proteins result 
in varying extents of entanglement in the protein hydrogels. As 
shown in Figure 2, hydrogels prepared via different chemistries 
all show an entanglement plateau modulus in the low-frequency 
regime. The chain molar mass distributions extracted from gels 
in Figure 1 do not fit well to the Jacobson-Stockmayer distribu- 
tion, likely because of challenges in controlling stoichiometry of 
reagents in small samples and because of difficulty in quantita- 
tively separating the large high molar mass band, which contains 
a mix of proteins with a molar mass of ca. 380kg/mol and 
above (6-mers and above). Therefore, the chain extension and 
the entanglement effect is quantitatively analyzed by the entangle- 
ment plateau moduli in the low frequency regime. The apparent 
entanglement molecular weight Mg can be estimated as (Larson 
et al, 2003) 

4 pwRT 
Me = -— — 

5 Ge 

with p being the protein density, cp the protein volume fraction 
in the gel, R the gas constant, T the absolute temperature and Ge 
the entanglement plateau modulus. Although the entanglement 
plateau modulus should not depend on the molecular weight of 
the protein polymers in the high molecular weight limit, proteins 
formed via macromolecular polycondensation show very broad 
molecular weight distributions with a certain fraction of the pro- 
tein chains below the critical entanglement length. These "small" 
proteins act as macromolecular diluents and lower the effective 
concentration of the entangled species, consistent with the con- 
centrations investigated being in the transition from the sticky 
Rouse phase to the sticky reptation regime (Tang et al., 2014). 
Therefore, increasing chain molecular weight, particularly the low 
molar mass tail, can increase Me- Consequently, differences in 
the degree of chain extension will yield differences in the plateau 
modulus of gels. 

Compared to disulfide bridging investigated previously (Tang 
et al., 2014), the three chemistries examined here show lower 
degrees of chain extension. Ge and Me are compared in 
Figures 3A,B- For 20% (w/v) hydrogels, disulfide coupling leads 
to a Ge of around 8800 Pa, approximately 2-8 times larger than 
the Ge's from other chemistries. As a result. Me is only 55 kg/mol 
in the disulfide coupling, even smaller than the molecular weight 
of a monomeric protein. The reasons for the lower entangle- 
ment molar masses in disulfide coupling are two-fold. First and 
most importantly, the stoichiometric imbalance between alkenes 
and thiols, due to inevitable experimental errors, set a practi- 
cal limit of the functional group conversion in the A-A -|- B-B 
type macromolecular polycondensation that becomes more acute 
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FIGURE 2 I Comparison of rheological frequency sweeps of chain The plateau modulus in the low frequency regime down to 0.001 rad/s is the 
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FIGURE 3 I Comparison of (A) entanglement plateau moduli, (B) the 
apparent entanglement molecular weights, and (C) the average 
entanglement densities. SS: disulfide, data from Tang et al. (2014); TM1, 



thiol-maleimide with la; TM2, thiol-maleimide with lb; TBM, thiol- 
dibromomaleimide; TE, thiol-ene. Error bars represent a 95% confidence 
interval. 



for chemistry performed on small samples. On the contrary, 
the disulfide bridging chemistry has no theoretical limit in the 
degree of chain extension due to stoichiometry since it is an 
A-A type polycondensation. Second, there are various side reac- 
tions in thiol-maleimide, thiol-bromomaleimide and thiol-ene 
coupling, which limit their ability to reach full conversion. As 
mentioned previously, the primary concern is that the occur- 
rence of side reactions may reduce the avaUability of the reactive 
functional groups (e.g., maleimide hydrolysis) or cause uncon- 
trolled chain coupling (e.g., amine -maleimide coupling). These 
side reactions further exacerbate challenges in controlling the 
stoichiometry of A and B reactive groups in the macromolec- 
ular polycondensation. Ge not only provides clues to analyze 
the extent of entanglement, but also provides information on 
the effect of side reactions. The difference in Ge might indi- 
cate that fewer side reactions occur in the thiol-ene coupling 
in the three chemistries examined in the study. Interestingly, 
although coupling with 2 only reaches moderate chain extension 
and side reactions might have the most detrimental effect, Ge 



from thiol-bromomaleimde coupling is larger those from thiol- 
maleimide chemistries. It is hypothesized that the supramolecular 
association between pCD and amino acid residues forms weak 
and dynamic crosslink junctions (with equilibrium constant log^T 
around 2-3) (Rekharsky and Inoue, 1998), thus increasing the 
modulus in the low-frequency regime. 

Creep experiments also provides information to compare the 
entanglement effects in gels prepared by different chemistries. 
Here, entanglement density E is used to quantify the entangle- 
ment effect, a quantity defined as the number of entanglements 
per molecule, namely, E = M/Me (Graessley, 2008). Strictly 
speaking, the definition above is used to describe the extent of 
entanglement of monodisperse polymers. Here the concept is 
borrowed to provide an estimate of the number-averaged entan- 
glement density. In the low E limit (when E < 6), the average 
entanglement density in the chain extended protein mixtures can 
be calculated as 

^avg — 2 . 5/^ Ge 
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where /° is the recoverable compHance, calculated as the extrap- 
olated intercept from a linear fit in steady state of the J-t data 
(Supplementary Figure 3); Ge is the entanglement plateau modu- 
lus, as defined by Graessley (2008) differently than the definition 
by Larson et al. shown above by a factor of 4/5. As shown in 
Figure 3C, the low E values confirm the previous assumption 
that the entanglement density is not large. It is also found that 
the average entanglement densities of the chain-extended pro- 
teins gives the same trend as the entanglement plateau moduli 
in gels, namely, disulfide > thiol-ene > thiol-bromomaleimide > 
thiol-maleimide. For entangled gels prepared via thiol-maleimide 
coupling, Eavg is only approximately unity. On the contrary, E„yg 
is larger than 3 for gels prepared via disulfide bridging. 

Because dibromomaleimide can be readily functionalized, 
dibromomaleimide coupling provides a facile method of incor- 
porating biological niches into the hydrogel network and embed- 
ding additional stimuli triggers to manipulate network structure. 
To illustrate this concept, dibromomaleimide is first derivatized 
with alkyne, and is further functionalized with pCD through 



copper-catalyzed azide-alkyne cycloaddition (CuAAC) to obtain 
2 (Scheme 2). pCD is chosen because its internal cavity can host 
many small guest molecules, such as drugs, via hydrogen bond- 
ing and hydrophobic interactions (Davis and Brewster, 2004). 
The most widely used complexation pair is pCD-adamantane 
with an association constant Ka of about 5 x 10* M^^ (Harries 
et al, 2005). This supramolecular association may offer modifica- 
tion sites to attach functionalities to the protein backbone. Here, 
the pCD-adamantane host-guest interaction is used to demon- 
strate the ability to add thermoreversible association to hydrogels. 
Adamantane end-capped monodisperse PNIPAM is synthe- 
sized by reversible addition-fragmentation chain transfer (RAFT) 
polymerization using trithiocarbonate 5 as the RAFT agent 
(Supplementary Figure 4), and the polymer's endgroup structure 
is confirmed by MALDI-TOF and NMR (Supplementary Figures 
5, 6, respectively). Upon complexation, the hybrid structure is 
hypothesized to adopt a graft-like structure: the chain-extended 
protein serves as backbone while the thermoresponsive polymers 
PNIPAM graft as side chains (Figure 4). Under the experimental 
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Scheme 2 | Synthesis of dibromomaleimide functionalized pCD (A) and adamantane functionalized trithiocarbonate RAFT agent (B). 
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conditions, approximately 91.8% of the PNIPAM can be attached 
to the proteins (see calculation details in Supplementary Materials 
and Supplementary Figure 7). The high-frequency plateau mod- 
ulus (G' at lOOrad/s) of the hybrid gel is lOkPa larger than the 
unmodified gel, and it increases moderately with increasing tem- 
perature, with a peak value at around 28° C. The entanglement 
modulus also increases from 2030 to 3670 Pa. The enhancement 
in the both moduli might originate from changes of protein 
architecture, and the formation of nanostructure within the gel 
due to segregation between the PNIPAM domain and the pro- 
tein domain. Such microphase separation may even exist at low 
temperature, as enhancement in moduli is observed as low as 
5°C. This observation is consistent with the previous finding that 
PNIPAM segregation takes place even when the polymers are 
highly solvated (Glassman and Olsen, 2013). 

INFLUENCE OF BRANCHED ARCHITECTURE ON THE MECHANICAL 
PROPERTIES 

Branched polymers have attracted much interest for fundamen- 
tal studies and industrial applications due to their distinct flow 
behaviors compared to their linear analogs (Dealy and Wissbrun, 
1990; Janzen and Colby, 1999; Wood-Adams and Costeux, 2001; 
Lohse et al, 2002; Graessley, 2008; Wang et al, 2010b). However, 
few reports characterize the synthesis or properties of branched 
proteins, using either chemical methods or cellular machinery 
(Zhang et al, 2013). Here, branched proteins are synthesized with 
the use of multifunctional crosslinkers in the chain extension 
reaction, demonstrating an important means to control the 
mechanical properties of protein gels by modulating the chain 
structure. The effect of branching is examined by varying the 



fraction of the trifunctional crosslinker (2,4,6-triallyloxy-l,3,5- 
triazine, 6) from 0 to 100% in the total alkenes (note that the 
ratio of total alkenes to thiols is kept 1:1). The hydrogels are pre- 
pared by thiol-ene click chemistry via thermal initiation. Similar 
to chain extension, all branching reactions are found to reach 
high endgroup conversion, yielding significant fractions of high 
molar mass proteins from SDS-PAGE (Figures). Significant dif- 
ferences in molecular weight distribution of proteins among the 
6 reactions examined cannot be measured by gel electrophore- 
sis. An exception occurs at 100% triene, where a larger fraction 
of low molar mass protein is present, suggesting that the end- 
group conversion might be lower than those in other cases. All 
protein hydrogels are able to be dissolved in 6M urea buffer 
post reaction, confirming the hypothesis that addition of 6 only 
results in branching and not crosslinking, even with 100% triene 
crosslinker. The gel point under each reaction condition can be 
calculated using the Carothers equation or the Flory-Stockmayer 
theory (see calculation details in Supplementary Materials and 
Supplementary Figure 8), which shows great dependence on the 
triene composition. The solubility of gels in urea suggests that the 
conversion of the endgroups is below the gel point in all cases. 

The mechanical properties of hydrogels are influenced by the 
branched structure of the proteins. To better compare the entire 
relaxation spectra with different fractions of 6, rheology data 
from creep experiments are converted to the dynamic compliance 
/ and /" using a Fourier transform (Ferry, 1980), 



/ (CO) =J°-w 
J" (w) = l/w)7o + CO 
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FIGURE 4 I Structure and mechanics of hybrid protein hydrogel. (A) Illustration of the brush-like structure of PNIPAM-extended-P4 chimeric molecule 
formed via pCD-adamantane complexation; (B) Temperature sweep at 1 °C/min heating/cooling rate; (C) Frequency sweep spectra. 



Frontiers in Chemistry | Chemical Biology 



May 2014 | Volume 2 1 Article 23 | 8 



Tang and Olsen 



Controlling entanglement in protein hydrogels 




B 0.25 



i 0.20 



0.00 El 



n — I — I — r-j — I — I — I — I — |— 1 — I — I — I — I — I — I — I — r 
0% 

5% 
10% 

20% A 
-C^ 35% 
100% 




1 2 3 4 5 6 

Degree of chain extension/branching 

FIGURE 5 I Molecular weight distribution of proteins in the branching 
reaction. (A) SDS-PAGE: The molar fractions of triene 6 in lanes 1-6 is 0, 5, 
10, 20, 35, and 100%, respectively; (B) Weight fraction extracted by 
densitometry. 



where /° is the recoverable compliance, co the angular frequency, 
)7o the zero-shear-rate viscosity, t the experiment time. /° and rjo 
are the intercept and the inverse of the slope, respectively, from a 
linear fit in steady-state phase of creep. The dynamic compliances 
are further translated into dynamic moduli G' and G", 

G' =//(;'' +/"') 

G"= 77 (/'+/"') 

The combined master curves are shown in Figure 6, which can be 
divided into three different regimes according to the relaxations 
present in the gels (from fast to slow): the coiled-coil regime, the 
entanglement regime and the terminal regime. With addition of a 
small amount of 6 (5 and 10%), both the high-frequency mod- 
ulus (G'oq) and the entanglement modulus (G^) decrease. This 
suggests that adding a small amount of crosslinker only causes 
short chain branching and small drops in the backbone length of 
the chain-extended proteins, since this would result in a decrease 
in the total number of entanglements. However, at 20% triene, 
both G'^ and Ge reach their peak values, where the number den- 
sity of the branch points and the branch length reach an optimal 
combination. In this case, the backbone length of the protein 
is likely to be comparable to the linear case, as the molecular 
weight distributions from SDS-PAGE of these two cases are not 
distinguishable, but the branched structure provides extra topo- 
logical interactions. The introduced branches retard the reptation 
of the entire molecule, as it is only after arm retraction that 
the reptation of the chain-extended backbone is allowed to hap- 
pen (Mcleish, 2002). Although the coiled-coil relaxation is only 
slightly affected by the molecule topology (see Supplementary 
Figures 9, 10), the reptation regime broadens in the frequency 
spectrum and extends toward the low frequency regime. This 



ro 10 
Q. 



10 r, 



10 



10 r 



_| 1 1 M MM| 1 


1 M Mll| 1 1 M nil| 1 1 M MM| 1 

entanglement : 
regime ^^g\ 


1 I 1 MM| 1 1 M llli| 1 1 1 1 Mil|. 


E terminal 


coiled-coil 1 
regime ; 


-_ regime^^ 








-•- 0% 
-•- 5% 
-•- 10% 
-•- 20% 
-•- 35% 
-•- 100% 


^ 1 1 


ii 1 ' 1 1 


Tl III 



10 10" 10 10" 10" 10 

Angular Frequency [rad/s] 



10 



10 




0% 5% 10% 20% 35% 100% 



Triene Composition 



FIGURE 6 1 Comparison of the mechanical properties of 
branched hydrogels at different triene compositions. (A) Master 
curves at 25°C from oscillatory shear and creep experiments; 
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the entanglement moduli, determined from the converted 
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is consistent with the observation of randomly branched poly- 
mers in the melt state (Mcleish, 2002). Further increase of the 
triene composition leads to decreases in both G'^ and Gg- This 
is hypothesized to originate from the decreased engroup conver- 
sion and the formation of hyperbranched molecular structures. 
The densely packed branches impose permanent topological bar- 
riers for different chains to entangle, which leads to much lower 
solution viscosities (Voit and Lederer, 2009). Recently, the vis- 
coelasticity and the dynamic relaxation of the synthetic hydrogels 
is found to be essential to mimic the complex biological tissues 
(Mckinnon et al., 2013). The presented branching reaction in this 
study establishes a facile approach of tuning the network mechan- 
ics over the entire frequency window, especially in the long time 
relaxation, which can be useful in mimicking the properties of 
natural tissues. 

CONCLUSION 

Under appropriate conditions, thiol-maleimide, thiol-bromo- 
maleimide and thiol-ene coupling can all result in significant 
chain extension. The differences in chemistries' reactivity and 
side reactions may cause variations in the molecular weight and 
its distribution of chain-extended proteins and the topological 
entanglement effect in gels. While thiol-maleimide conjugation 
is the most common and the easiest to implement, thiol-ene click 
chemistry can achieve high endgroup conversion fairly rapidly. 
Thiol-bromomaleimide shows its potential in reversibly modi- 
fying proteins, and using functionalized dibromomaleimide as a 
chain extender allows further control of hydrogels' properties by 
incorporating side chain functionalities into the protein architec- 
ture. Here, thermoresponsive changes in mechanical properties 
of gels are demonstrated with PNIPAM grafts. Lastly, branched 
proteins are prepared in A2 -|- B2/B3 mixed type reactions. The 
entanglement plateau modulus is increased when branches are 
long enough to enhance topological constraints yet the added 
branches do not sacrifice the backbone length. In conclusion, 
the structure of engineered proteins and their assembly behav- 
iors can be easily modified with the use of different chemistries, 
which presents a rich toolbox to tailor the structure-properties of 
protein materials in various applications. 
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